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ABSTRACT 



Aims. We present an analytic approximate seismic inversion scheme for damped transverse coronal loop oscillations based on the thin 
tube and thin boundary approximation for computing the period and the damping time. 

Methods. Asymptotic expressions for the period and damping rate are used to illustrate the process of seismological inversion in a 
simple and easy to follow manner. The inversion procedure is formulated in terms of two simple functions, which are given by simple 
closed expressions. 

Results. The analytic seismic inversion shows that an infinite amount of 1 -dimensional equilibrium models can reproduce the observed 
periods and damping times. It predicts a specific range of allowable values for the Alfven travel time and lower bounds for the density 
contrast and the inhomogeneity length scale. When the results of the present analytic seismic inversion are compared with those 
of a previous numerical inversion, excellent agreement is found up to the point that the analytic seismic inversion emerges as a 
tool for validating results of numerical inversions. Actually it helped us to identify and correct inaccuracies in a previous numerical 
investigation. 
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1. Introduction 

References to coro nal seismology date back to the 1 980s 
([Roberts et al.1 11984 and even the 1970s (lUchidal [197 0^, but 
coronal seismology remained largely a theoretical concept. This 
situation changed drastically when observations from space ob- 
servatories showed that MHD waves are ubiquitous in the so- 
lar atmosphere. Opinions might differ, but we are inclined to 
identify the detect ion of damped transvers e co ronal loop oscilla- 
tions i n 1999 bv Aschwanden et alld 19991) and lNakariakov et al] 
(1 19991) in observations made with the EUV telescope on board 
of the Transition Region and Coronal Explorer (TRACE) as the 
real start of coronal seismology. Since the n the detection of thes e 
oscillations has been confirmed (see e.g. ISchriiver et al.ll2002l) . 
The TRACE oscillations have periods T of the order of ^ 2 - 10 
minutes and comparatively short damping times of the order of 
Td - 3 - 20 minutes. There is general consensus that the TRACE 
oscillations are fast standing kink mode oscillations. In addi- 
tion, damped oscillations observed in hot coronal loops by the 
SUMER instrument on board SO HO have been interprete d as 
standing slow mode oscillations dWang et al.ll2002i l2003h and 
the measured period has b een used to det ermine the magnetic 
field strength in the loop ( Wang et al.ll2007i) . 

Theory shows that when the fast magneto-sonic kink MHD 
waves have their frequencies in the Alfven continuum, they cou- 
ple to local resonant Alfven waves (W right & Rickard 1995; 
iTirrv & Goossens|[l996l: iGoossens et al. i r2006t lGoossensll200i 
and get transformed into quasi-modes that are damped by reso- 



nant absorption dTirrv & Goossenslll996l) . This is exactly what 
happens for a radially stratified cylindrical coronal loop model 
(Goossens et al. 2002; Ruderman & Roberts 2002). Coupling of 
the fast magneto-sonic MHD waves to local Alfven waves is a 
natural phenomenon in non-uniform coronal loops. This reso- 
nant coupling produces a quasi-mode which in a static equilib- 
rium is damped by resonant absorption. The damping is inde- 
pendent of the dissipative coefficients. 

This finding does not mean that other damping mechanisms 
are ruled out. Resonant absorption might not be the only cause of 
the observed damping but it is definitely operational. Strong sup- 
port for the robust character of quasi-modes and their resonant 
damping comes from a recent investigation by iTerradas et al.l 
(2008) on MHD waves in multi-stranded coronal loops. An im- 
portant finding of this investigation is that resonantly damped 
quasi-modes live on complicated multi-stranded coronal loops. 
They do not need nice cylindrical magnetic surfaces as might be 
concluded from studies on simplified 1 -dimensional equilibrium 
models. The message is that 1 -dimensional models are a great 
help to reduce the mathematical complexity but still contain the 
essential physics of resonantly damped quasi-modes. 

iRuderman & RobertsI (l2002h were the first to suggest that the 
observed rapid damping of the transverse oscillations of coronal 
loops could be explained by resonant absorption. In the con- 
text of the heating of solar plasmas HoUweg &JYang (19881) 
have predicted that oscillations in c oronal loops are to un dergo 
rapid damping. In the same context iGoossens et al.l (Il992i) de- 
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rived analytical expressions for the frequency and the damping 
rate of quasi-modes in static and stationary equilibrium models. 
[Ruderman & Roberts ( 2002) focused on proving the principle 
of resonant absorption as damping mechanism for the transverse 
oscillations i n coronal loops and co nsidered one specific numeri- 
cal example. lGoossens et al.l (|2002|) looked at the damping times 
of 11 loop oscillation events and basically confirmed that reso- 
na nt absorption can explain the observed damping as suggested 
bv lRuderman & RobertsI (l2002h . 

The analytic expressions derived bv iHollweg & Yand 

([1988): Goossens et all (1 19921) : iRuderman & Roberts! (|2002|) . 
and iGoossens et al.i (12001 for the damping rate are asymptotic 
in the sense that they are derived in the assumption that the 
non-uniform layer is thin. This is the so-called thin boundary 
approximation, in what follows we shall refer to it as the TB- 
approximation. 

The seismologica l studi es on transverse oscillations 
are bv .Nakariakovl (120001) 
' (120021) : lAschwand en et al 

iNakariakovi ( l2000h and 
(1200 Ih used the observed periods 



so far 

dlOOl : 



iGoossens et al 



(2007). 



and 'Arreg ui et al.l 
[Nakariakov & OfmarJ 
and theoretical estimates of the periods 
wavelength or thin tube approximation 



iNak ariakov & Ofmanl 
(120031) . 



based on the long 
(TT-approximation) 
for a uniform coronal loop model to derive estimates for the 
magnetic field strength. The weak link in their analysis is the 
uncertainties on the density. Goossens et al. (2002) used the 
observed damping rates and theoretical values of the damping 
rates based on the TB-approximation to derive estimates for 
the radial inhomogeneity length -scale. Again, the weak l ink is 
the uncertainties on the density. iAschwanden et al] (l2003l) used 
the observed damping rates and the damping rates computed 
by IVan Doorsselaere et al.l (12004 . outside the TB -regime, to 
determine the density contrast. The first study that used the 
observational information on both periods and damping times 
in the context of res onant damping in a consistent manner is by 
lArregui et al] (|2007|) . 

The important finding of lArregui et al.l (l2007b was that an 
infinite amount of 1 -dimensional cylindrical equilibrium mod- 
els can reproduce the observed period and damping rate with 
the internal Alfven transit time or conversely the int ernal Alfven 
veloci ty confined to a short range. The study by .Arregui et alj 
( l2007h is fully numerical and apparently part of the physics has 
remained not well understood. The objective of the present paper 
is to use asymptotic expressions for the period and damping rate 
to illustrate in a simple and easy to follow manner the process 
of seismological inversion. The asymptotic expressions are the 
TT-approximation for the period and the TB-approximation for 
the damping rate. When both approximations are used simulta- 
neously we shall refer to it as the TTTB-approximation. We are 
well aware of the fact that in case of strong damping this approx- 
imation might give inaccurate results. However, our primary in- 
tention is to understand the process of seismic inversion and as 
we shall see the asymptotic expressions turn out to be accurate 
far beyond their theoretical range of validity. 

2. Asymptotic analytic expressions for period 

The analytical expression that we shall use for the period is ob- 
tained by (i) adopting the TT-approximation for MHD waves and 
by (ii) modelling the coronal loop as a uniform cylinder with a 
straight magnetic field along the z-axis. The TT-approximation 
means that the period is independent of the radius and that ef- 
fects due to non-zero radius are absent as far as the period is 
concerned. The choice of a uniform equilibrium model means 



that effects due to stratification are absent. The coronal loop is 
modelled as a cylindrical plasma with constant density pi embed- 
ded in an external plasma with constant density p^. The coronal 
loop is basically a density enhancement with pi > p^. The mag- 
netic field is constant and has the same strength both inside and 
outside the loop. 

Our starting point is the well-known expression for the 
square of the frequency of the kink mode in a uniform cylin- 
der with a straight ma gnetic field along the z-axis (see e.g. 
lEdwin & Robertsll 19831) 



Pi +Pe 



(1) 



The subscripts i and e refer to quantities respectively in the 
coronal loop and in the external plasma surrounding the loop. 
a>A = kzVA is the local Alfven frequency, with the parallel 
wavenumber, Va = B/ VO"P) the local Alfven velocity, and B 
the magnetic field strength. Hence we can rewrite Eq. ([T]i as 



B 



(2) 



with ( - pi/pf. > 1 the density contrast. Now we note that for 
the observed transverse oscillations, with a wavelength double 
the length L of the loop, - njL and we convert frequenc ies to 
periods and rewrite Eq. (|2| as (see e.g. lArregui et alj[2008l) 



T = TA,, V2 AiO- 



(3) 



Here TA,i = LIVa^ is the internal Alfven travel time and the 
function A (f) is defined as 



A{0 



^+1 



1/2 



(4) 



Equation (O is our first key equation. Let us recall that 
Eq. (O has been obtained by use of the TT-approximation, hence 
effects from non-zero radius and stratification are absent. This 
equation expresses the period T, which is an observable quan- 
tity, in terms of the Alfven travel time tax and the density con- 
trast ( which are two quantities that we aim to determine from 
seismic inversion. If we have observed values of the period T 
and we convince ourselves that Eq. (O is a good first analytical 
approximation of the period T then we can invert it for either TA,i 
or (. Actually, we shall do both. Let us first solve Eq. ^ for ta,!- 
Since we prefer to use dimensionless quantities we introduce y 
as 



y = 



TA,i 



From Eq. (|3]l we obtain 
T 



1/2 



(5) 



(6) 



Since we have not any information on ^, it might appear that 
Eq. (|6]l is not very helpful. However, closer inspection reveals 
that it contains very valuable information. For a given observed 
period T, Eq. (|6|l is a parametric representation of y (or equiva- 
lently of TA,i) in terms of (. In order to stress this point we define 
the function F\ by use of the right hand member of Eq. ^ as 



[1, oo[ ^ R, ( 



V2U+I 



1/2 



(7) 
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For later comparison in Sect. |4] we note that dFi /d^ > and 
d^Fi /d(^ < 0. Hence y = Fi is a strictly increasing and con- 
cave function of We also note that 

This means that 



y < 



1 



T 



(8) 



According to inequality ((SJ the Alfven travel time is constrained 
to a narrow range; whatever the density contrast is, the Alfven 
travel time is in between 0.5 T and 0.707 T . If we accept that 
the density contrast ^ is not smaller than say 1.5 then 0.5 T is 
replaced with 0.548 T narrowing further down the range for ta.i. 

If we are able to determine the length L of the loop, then we 
can extract Va^ and find 



Let us recapitulate what seismic information we have de- 
duced from the observed value of the period. The quantities ( 
and y are in the following intervals 



and are related to one another by 



(14) 



(15) 



The functions F\ and G\ are defined by expressions (flJl i. Of 
course, G\ is the inverse function of F\. G\ - ^ ' 



' and con- 



versely G[ - Fi. When the period T is known from observa- 
tions, then there are infinitely many pairs (f,y) that reproduce 
the observed period. We can let ^ vary over the interval [1, oo[ 
and compute for each value of ( the corresponding value of 
y - F\ (f ), or conversely, let y vary over the interval [1/2, 1 / V2[ 
and compute for each value of y the corresponding value of 
^ = Gi(y). 



Va,; = V2 - AiO, 

r L L 
V2 - < Va,! < 2-. 



(9) 



Hence, here is the na rrow range for the local Alfven velocity that 
lArregui et al.l (l2007l) found in their numerical inversion. If the 
density contrast is not smaller than 1.5, 2 LjT is replaced with 
1.825 LIT further reducing the available range for Va^ to less 
than a 25% relative margin compared to its possible maximal 
value. 

With the help of the first Une of inequality (O we can refine 
the definition of F\ and replace (Q with 



[1, 



Let us now solve Eq. (|3]l for ( and find 
l-2y 



1/2 



,2 ■ 



.(10) 



(11) 



Equation ( fTTT i is the twin of Eq. For a given observed period 
T, Eq. ( fTTI ) is a parametric representation of ^ in terms of y (or 
equivalently in terms of taj). In order to stress this point we de- 
fine the function G\ by use of the right hand member of Eq. ( fTTT l 

as 



2 V2 



y 



Gxiy) 



2/ 



1-2/ 



(12) 



Since dGi /dy > and d^G i /dy^ > 0, it follows that ^ = G i (y) is 
a strictly increasing and convex function of y. Note now that 

Gi(l/2) = 1, lim Gi(y) = oo. 

v^l/ V2 

With this information on the function Gi we can refine its def- 
inition ( fT2] i. Combined with the definition ( fTOl l of F\ we obtain 
the following prescriptions for the functions F\ and G\ ; 



Fi : [1, <x,[^ [i, -^l ^ ^ FdO 
2 V2 

Gi : ^ [1, y ^ Gi(y) = 

2 V2 



V2U+I 

2y^ 



1/2 



1 -2y^ 



(13) 



3. Asymptotic analytic expressions for damping 
time 

In order for the kink MHD waves to be damped by resonant ab- 
sorption additional physics has to be introduced in the equilib- 
rium model. The required additional physics is non-uniformity 
of the local Alfven velocity. For a constant magnetic field 
this implies a non-uniform density. Asymptoti c expressions for 
the damping time have been deri ved bv e.g. |Hollweg & Yand 
( ll988h . lGoossens et al.l ( ll992l) and lRuderman & Robertsl ( l2002lK 
These asymptotic expressions are derived in the approximation 
that the non-uniform layer is thin. The true density discontinuity 
is replaced by a continuous variation in density. Jump condi- 
tions are used to connect the solution over the ideal singularity 
and to avoid solving the non-ideal MHD wave equations. The 
jump condition for the id eal Alfven s ingularity was introduced 
on an intuitive manner bv lHollweg & Yang ( 198 8) and put on a 
firm rn athe matical basis by Sakurai et a l. ( 1991). lGoossens et al.l 
(Il995h and|Goossens & Ruderman' ('1995) for the driven prob- 
lem, and bylTirrv & Goossens (1996) for the eigenvalue prob- 
lem. The result of this asymptotic analysis is 



II - f ^ Pi """Pe 
T I Pi- Pi 



(16) 



Here T and are the period and the damping time, pi is the in- 
ternal density on the axis of the loop and hence in the interval 
[0, R - 5]; pe is the constant external density where external 
refers to[R+ i^, oo[, / is the thickness of the non-uniform layer; 
II R - corresponds to a uniform loop and a discontinuous vari- 
ation in the density at the radius R of the loop. A fully non- 
uniform loop has IIR - 2 but cannot be studied by an approxi- 
mate TB expression for the damping time. The numerical factor 
F depends on the variation in density across the non-uniform 
layer. For a linear variation as used b y iGoossens et alj (Il992h 
F — 4/7r^. lRuderman & RobertsI (l2002 f) used a sinusoidal varia- 
tion in density across the non-uniform layer and found F = lln. 
In what follows we shall adopt F - 2 In. A sinusoidal varia- 
tion in density in the non-uniform transitional layer is proba- 
bly closer to reality than a linear variation. In addition we shall 
compare our analytic results with numerical results for fully 
no n-uniform 1 -dimensional mod els of cor onal loops o btained 
by IVan Doorsselaere et all ( 12004) and by Arregui et alj (l2007h 
by the use of the eigenvalue code LEDA originally designed by 
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I Van der LindenI (Il991h . In these studies a sinusoidal variation in 
density was considered. Note also that in Eq. ( fT6l ) the period is 
computed using the TT-approximation. 
Rewrite Eq. ( fTSI l in terms of ^ to find 



T 



2^+1 1 
n(-l JJr' 



(17) 



Equation ( fTTj l is the second key equation of the present inves- 
tigation. Let us now look at this equation from a seismic point of 
view. If we have observed values of the period T and the damp- 
ing time Td and we convince ourselves that Eq. ( [TT] ) is a good 
first analytical approximation of the damping time divided by 
period then we can invert Eq. (fTTT i for either ( or l/R. Actually 
we shall do both. Let us first solve Eq. ( [TT] ) for ^ and find 



— — + 1 

2R T ^ 

I TTTd 

JR T 



1 



(18) 



From a seismic point of view, period and damping time can be 
considered as known from observations and it makes sense to 
denote tit^IT as a constant C. In addition it is convenient to 
abbreviate l/(2R) as z- Hence 



I 

-R' ^ 



Equation (fTSl l can be rewritten as 



Cz + 1 
Cz-V 



(19) 



(20) 



Here are several observations to be made. First, since ^ > it 
follows from Eq. (|20] | that for a given ratio of T^jT there is a 
lower bound for the inhomogeneity length scale 



_ I r _ 1 _ 

Z — ■—— > — — — Zmin- 

2R TTTd C 



(21) 



Second, realize that Eq. ( |20] i is a parametric representation of ( 
in terms of z = 1/{2R). In order to make this point very explicit, 
we introduce the function G2 defined as 



G2 :] 



1] ^ R, z ^ G2(z) = 



Cz+l 
Cz-l 



(22) 



It is easy to show that dGi/dz < and d^G2/dz^ > 0. Hence 
^ — Giiz) is a strictly decreasing and convex function of z- It 
attains its absolute minimum for z = ll(2R) - 1 . This minimal 
value is 



on the periods. If we also use the information on the damping 
rate then we obtain 



y 



With the help of the information on the bounds for y and ^ 
we can refine the definitions ( fT3] l for F\ and G\ respectively to 
their final versions 



c + 1 Alc + i 

C-1 ^ ^2\ C 



1/2 



1 

V5 



1/2 



Gi 



1/C^ 
^2 I C 



1/2 



y 



Giiy) = 



V2 

2/ 



z - 



1-2/- 

Let us now solve Eq. ( fTTl i for z-lj (2R) and find 

/ _ 1 ^ + 1 
2R ~ C^-V 



(26) 



(27) 



Equation (|27] | is a parametric representation of z = 1/{2R) in 
terms of As before we make this point explicit by introducing 
the function F2 defined as 



(28) 



It is easy to show that dF2/d( < and d^F2/d(^ > 0, hence 
z - F2iO is a strictly decreasing and convex function of ^. In 
addition 

i^2(^mm)=l, limF2(0=^- 

With this information on the function F2 we can refine its def- 
inition (|28] |. Combined with the definition (l24l l for the function 
G2 we obtain 



F2 : [ 



,Ch- 1 

c- r 



4^] 



1], f ^ F2(0 



G2 : ]-, 1] 



[ , °°l 



G2(Z) 



1 ^+1 

c^- r 
_ cz-H 1 

" Cz- 1' 



(29) 



Here F2 is the inverse function of G2', F2 = Gj' and conversely 
G2 is the inverse function of F2- 



^min - G2(l) 



c-i' 



(23) 



Conversely ( attains its maximal value of infinity in the limit 
z — > 1/C. With this information on the function G2 we can refine 
its definition (l22l) as follows 



G2 : 



1] 



X- 1 



G2(Z) 



Cz+1 
Cz-l' 



(24) 



With the minimal value ( |23] l for ^ we can slightly improve on 
the lower bound 1 /2 for y and T/2 for ta.i given by inequality (|8j 
and on the upper bound 2L/T for VA,i given by inequality (|9]). 
These bounds were found with the sole use of the information 



4. Analytical seismology and comparison with 
numerical inversion 

Let us recapitulate the key results of the previous section. The 
two quantities that we assume to be known from observations are 
the period T and the damping time Td. Analytical theory based 
on the TTTB-approximation gives us two equations, namely 
Eqs. (O and ( [Tt] ) that express the period T and the damping time 
Td in terms of the density contrast the Alfven transit time (nor- 
malised to the period) y - TA,i/T and the inhomogeneity length 
scale (normalised to the radius of the loop) z - ll(2R). These 
three quantities ^, y and z are the seismic quantities in the sense 
that they are the quantities that we aim to determine with the 
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Loop 


R 


L 


R/L 


T 


Td 




TA.iA 


TA.iN 


Va.IN - L/TajN 




(10'^ m) 


(10'* m) 


(10-2) 


(s) 


(s) 




(s) 


(s) 


(km s-') 


1 


3.60 


1.68 


2.13 


261 


870 


0.30 


137 (143)~185 


145-177 


947-1161 


2 


3.35 


0.72 


4.65 


265 


300 


0.88 


150-187 


163-182 


396-443 


3 


4.15 


1.74 


2.37 


316 


500 


0.63 


173-224 


189-217 


801-922 


4 


3.95 


2.04 


1.94 


277 


400 


0.69 


153-196 


168-189 


1079-1211 


5 


3.65 


1.62 


2.26 


272 


849 


0.32 


143 (150)-192 


151-187 


869-1073 


6 


8.40 


3.90 


2.17 


522 


1200 


0.44 


279 (286)-369 


304-359 


1088-1284 


7 


3.50 


2.58 


1.37 


435 


600 


0.73 


241-308 


267-299 


862-967 


8 


3.15 


1.66 


1.88 


143 


200 


0.72 


79-101 


90-98 


1693-1853 


9 


4.60 


4.06 


1.15 


423 


800 


0.53 


229 (232)-299 


247-291 


1394-1643 


10 


3.45 


1.92 


1.78 


185 


200 


0.93 


105-131 


117-126 


1526-1643 


11 


7.90 


1.46 


5.38 


390 


400 


0.98 


227-280 


245-270 


541-595 



use of observed values of the period T and the damping time Td. 
Since we have only two equations that relate the three unknown 
quantities to the two observed quantities there are an infinite 
number of solution s of the seism i c inve rsion problem, such as 
first pointed out bv lArregui et al.l (l2007h . The seismic variables 
are constrained to the following intervals 

z € 1] (30) 

and are related to one another by 

y = F,(0. ( = Gxiy), 

z = F2(0, ^ = G2(z). (31) 

The functions F\,G\,F2,G2 are defined by ( |26] l and ( |29] l. 

Of the four functions only two are independent since G\ is 
the inverse function of F\ and G2 is the inverse function of F2- 
The set ( [3TT l gives us the infinitely many solutions of the seis- 
mic inversion in parametric form. Each of the three unknowns 
can be used as parameter and the two remaining unknowns can 
be expressed in terms of that parameter For example choose 
^ as parameter. Let ^ take on all values in and compute the 
corresponding values of y and z by the use of y - Fi(() and 
z = F2{0- Or choose y as parameter. Let y take on all values in 
/,. and then compute the corresponding values of ( and z by the 
use of ^ = Giiy) and z = ^2(^1(3')). Finally, use z as parameter 
to define the solutions of the inversion problem. Let z take on all 
values in /, and then compute the corresponding values of y and 
^ by the use of ^ = G2(z) and y - Fi(G2(z)). 

As an illustrative example we re-anal yse loop oscillation 
event # 5 examined by -Arregui et alj (12007 ) in detail using their 
numerical seismic inversion scheme. The results of the investiga- 
tion of that loop event are shown on their Figure 3a. For this loop 
oscillation event T - 272 s and Td - 849 s. The radius R and 
the length L of this loop are estimated to be /? = 3.65 x 10'' m 
and L - 1 .62 x 10^ m respectively. The ratio of the radius to the 
length of the loop is R/L ^ 2 x 10"^. This small value is good 
news for the TT-approximation to the period with effects due to 
a non-zero radius on period being small. 

The constant C-{nTd)/T = 9.81. Hence ^min - 1.23, y^in - 
0.525, jniax = 0.707 and z^in = 0.102. The intervals for ^,y,z 
are now 

( e [1.23, oo[. 



y € = [0.525, 0.707 [, 

z e /, = ]0.102, 1]. (32) 

The corresponding interval for TA,i is 

143 s < TA,i < 192 s. 

If we limit the analysis to ^ > 1.5 then the lower bound is 150 s 
so that 

150 s < TA,i < 192 s. 

This c an be compared with the interval in Fig 3a of lArregui et al.l 
(I2007h which is (see also their Table 1) 

170 s < TA,i < 210 s. 

It is encouraging to note that the overall differences on the 
Alfven travel times found here are about 10% although the loop 
oscillation event is characterised by heavy damping. It is intrigu- 
ing that the interval as a whole is shifted by about 20 s to shorter 
Alfven travel times. This intriguing discrepancy has led us to 
calculate analytical results for the allowable range of the Alfven 
travel time for the remaining 10 l oops a nd to compare them with 
those obtained bv lArre^ui et alj (l2007h . It turned out that in all 
cases there are discrepancies and some of them were far bigger 
than 10%. This has motiva ted us to re-examine the numerical re- 
sults of the investigation of Arregui et alj (|2007[) . It turns out that 
the values given in Arregui et al. (2007]) are inaccurate. We have 
re-calculate d the numeric al values for the Alfven travel time 
intervals in lArregui et al.l (12007.) and present them in Table [T] 
When this issue is taken into account, the corrected numerical 
interval for the event under consideration is 

151 s < TA,i < 187 s, 

which remarkably agrees with the analytic interval. Table [T] 
shows that analytic and numerical Alfven travel time intervals 
agree very well for all events. Analytic lower bounds for TA,i are 
a little below numerical ones. Note that analytic lower bounds 
for ^ given in ( [30l l can be slightly lower or larger than the nu- 
merical f = 1.5. When lower, values in parentheses give the 
analytic TA,i for this density contrast. Upper bounds correspond 
to the limit TA,i(^ — > 00) = r/ V2, which is nearly approximated 
at the numerical ( ^ 20. We have also calc ulated Zmin = 1 /C for 
the events studied by 'Arregui et ajj (|2007[) and the results agree 
well with the corresponding values of (l/R)mm in their Table 1 . 

As an illustration that any of the three seismic quantities can 
be used as parameter we take z as the parameter We let z vary 
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Table 2. Analytic seismic inversion results for loop # 5. 



z = l/(2R) 


( 


y = tm/t 


0.105 


67.67 


0.702 


0.110 


26.32 


0.694 


0.120 


12.30 


0.680 


0.125 


9.85 


0.674 


0.150 


5.24 


0.648 


0.175 


3.79 


0.629 


0.200 


3.08 


0.617 


0.225 


2.66 


0.603 


0.250 


2.38 


0.593 


0.275 


2.18 


0.585 


0.300 


2.03 


0.579 


0.325 


1.91 


0.573 


0.350 


1.82 


0.568 


0.375 


1.75 


0.564 


0.400 


1.68 


0.560 


0.425 


1.63 


0.557 


0.450 


1.59 


0.554 


0.475 


1.55 


0.551 


0.500 


1.51 


0.549 


0.525 


1.48 


0.546 




and compute for each value of z the corresponding value of ^ by 
the use of the function G2 defined in ( |29] l. 



9.81Z+ 1 

9.8iz- r 



(33) 



Subsequently we compute the corresponding value of y by the 
use of the function Fi defined in i 



y 



V2U+1 



1/2 



(34) 



Of course, we can only compute ^ and y for discrete val- 
ues of z- The results of our computations are summarised on 
Table |2] and graphically represented on Fig. [T] Recall that the 
functions Fi,G\,F2,G2 (see eqs. 1261 and 1291 ) are monotoni- 
cally increasing (F\,G\) or monotonically decreasing (^2,^2) 
and have concave graphs {F\) and convex graphs (F2,GuG2) 
respectively. This implies f.e. that y = TA.ilT is a strictly in- 
creasing function of ^ and a strictly decreasing function of z and 
conversely that z is a strictly decreasing function of both ^ and 
TA.i- Inspection of the second order derivative of a given quan- 
tity with respect to one of the two remaining quantities can tell 
us that the graph is either concave or convex. For example the 
graphs of y and of z as function of ( are respectively concave 
and convex. The monotonic variation of ^, y, z and the concave or 
convex appearance of their graphs predicted by our analytic seis- 
mic inversion agrees exactly with the behaviour of the numerical 
inversion. Furthermore, Fig. [T]displays an amazing quantitative 
agreement between analytic and numerical inversion results. 

We have not re-ana lysed in detail l oop oscillation event #10 
that was examined by lArregui et alj (12007 ). The inversion for 
that loop event is shown on their Figure 3b. A striking property 
of the solutions is the non-monotonic behaviour of the seismic 
variables. This is clearly reflected in the pronounced minimum 
of TA,i as function of ^ and as function of IjR. The decreasing 
part of TA,i as function of ^ and the increasing part of taj as 
function of IIR cannot be recovered by the analytical seismic 
inversion scheme based on the TTTB-approximation. The ana- 
lytical TTTB-approximation predicts monotonic variation of the 
seismic variables and cannot approximate multi-valued solutions 
with two pairs of (^, IjR) corresponding to the same value of ta.i. 



Fig. 1. Analytic inversion (solid lines) and corrected numerical 
inversion (filled circles) in the (^, IjR, TA,i)-space for loop oscil- 
lation event #5 in Table [T] 



The fact that the analytical seismic inversion fails for this loop 
oscillation event is not disturbing since this event is characterised 
by extremely heavy damping with T/rd - 0.92 which we antic- 
ipated would fall out of the application range of the analytical 
scheme anyway. 



5. Conclusion 

In this paper we have presented an analytic approximate seis- 
mic inversion scheme based on the TTTB-approximation for 
computing the period and the damping time. In the TTTB- 
approximation the period is computed for a uniform loop model 
in the long wavelength or zero radius approximation. The damp- 
ing time is computed for relatively weak damping corresponding 
to thin non-uniform layers. The advantage of this analytical seis- 
mic inversion is that it is formulated with the aid of two functions 
Fi and F2 (and their inverse functions G\ - F^'^ and G2 - F^^) 
which are given by simple closed expressions. The practical im- 
plementation of the inversion scheme is stunningly simple. The 
calculations required to obtain solutions can even done with the 
use of a hand calculator. This analytical scheme seismic inver- 
sion clearly shows that the inversion problem has infinitely so- 
lutions in the (^, y, z)-space as first pointed out by lArregui et al.l 
(120071) . It also reveals that the allowable values of y (or Alfven 
travel time) are confined to a narrow range. When applied to a 
loop oscillation event with heavy damping as f.e. loop oscilla- 
tion event # 5 with T jr^ - 0.32 the analytic inversion scheme 
produces remarkably accurate results. Not only does it recover 
the overall appearance of the solution curve with the correspond- 
ing monotonic behaviour of the seismic variables. In addition, it 
recovers for a prescribed range of values of ( the corresponding 
values of y (or taj) and z (or IIR). 

The disadvantage of the scheme is that (i) it does not take 
into account the effects of non-zero radius and of radial stratifica- 
tion on the period; (ii) it is strictly speaking only valid for weak 
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damping corresponding to thin non-uniform layers. Corrections 
due to finite tube radius are of the order of the loop radius to 
length ratio squared. For the largest observed value of RjL in 
Table[T]the correction is of the order of 10"^, hence the thin tube 
approximation does not impose any practical restriction on the 
applicability of analytical results. As for the thin boundary ap- 
proximation, for cases with extremely heavy damping, the non- 
monotonic behaviour displayed by numerical solutions cannot 
be recovered. This is the price for an analytical scheme. All in 
all, the accuracy of the results obtained with this analytic in- 
version is amazing. The final agreement of the analytic seismic 
inversion with the numerical seismic inversion when the inaccu- 
racies of the numerical inversion are removed is excellent up to 
the point that the analytic seismic inversion emerges as a tool for 
validating results of numerical inversions. 
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